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PROCEDE DE FILTRAGE DE DONNEES SiSMIQUES. NOTAMMENT 

PAR ANALYSE KRIGEANTE 

DOMAINE TECHNIQUE GENERAL - RAPPELS SUR L'ANALYSE KRIGEANTE 

5 

La presente invention est relative au filtrage de donnees sismiques, 
notamment par analyse krigeante 

L'analyse krigeante permet la decomposition d'une fonction 
aleatoire a partir de sa fonction de covariance. 

10 Elle est notamment classiquement utilisee en geostatistique pour 

filtrer des donnees sismiques, notamment, mais non limitativement, pour 
caracteriser des r6servoirs. 

L'analyse Krigeante repose notamment sur I'hypothese qu'un 
phenomene mesure localement au moyen d'un echantillonnage regulier ou 

15 non peut s'analyser comme la somme lineaire de plusieurs phenomenes 
independants, le variogramme du phenomene global correspondant a la 
somme lineaire des variogrammes de chacun des phenomenes 
independants qui le composent. 

Classiquement, on decompose le variogramme correspondant aux 

20 donnees experimental mesurees en une somme de variogrammes 
modelises et on deduit des donnees experimental et des modelisations 
choisies pour les variogrammes elementaires de la decomposition, les 
fonctions elementaires qui decomposent la fonction aleatoire correspondant 
au phenomene global. 

25 II est ainsi possible d'extraire d'une cartographie de donnees 

sismiques du type de celles illustrees sur la figure 1 (par exemple, des 
donnees experimental brutes) d'une part le bruit blanc present dans ces 
donnees (figure 2a), d'autre part un bruit correspondant a des raies lineaires 
(figure 2b) et enfin des donnees filtrees, debarrassees de I'un et I'autre de 

30 ces deux bruits (figure 2c). 

Les calculs de determination par krigeage des valeurs des fonctions 
elementaires decomposant une fonction aleatoire globale sont en eux- 
memes classiquement connus par I'homme du metier. 
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On pourra par exemple a cet 6gard se referer aux articles et 
publications mentionn6es dans la bibliographie jointe a la fin de la presente 
description. 

De fagon tres generate, on determine la valeur d'une fonction 
5 el6mentaire intervenant dans la decomposition de la fonction aleatoire 
globale comme etant une combinaison lin£aire de valeurs exp6rimentales 
pour des points d'un voisinage immediat du point consider6, ces valeurs 
exp6rimentales 6tant affectees de coefficients de ponderation. 

En d'autres termes, si Ton considere qu'une fonction Z(x) se 
10 decompose en la somme de fonctions elementaires Y u (x), 
ce qui se traduit par : 



ou a est un indice muet designant les points consid6res autour du 
20 point x pour lequel on cherche a determiner la valeur estimee Y u *(x), Z x 
6tant la valeur au point x, N etant le nombre de ces points. 

On demontre que les coefficients de ponderation Xa verifient 
Tequation : 



15 



on estime la composante Y u (x) par : 



Y u *(x)= ± fa* 




^01 



25 



oCi I'indice 0 designe le point pour lequel on cherche a determiner 
une estimation, les valeurs C 0 " a C% N etant des valeurs de covariance 
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calculees a partir du modele u correspondant a la composante Y u (valeurs 
de la fonction de covariance pour les distances entre chaque point de 
donnee et le point a estimer), les valeurs Cy etant des valeurs de 
covariance calculees en fonction de la mod6lisation choisie pour le 
5 variogramme de la fonction & estimer (valeurs de la fonction de covariance 
pour les distances entre les points de donnees). 

Ces coefficients de pond§ration X x u se d£terminent done par simple 
inversion des matrices de covariance. 

10 PROBLEME(S) POSES PAR LET AT DE LA TECHNIQUE - PRESENTATION DE 
L'INVENTION 

Une des difficult^ des techniques d'analyse Krigeante connues a 
ce jour est qu'elles n6cessitent d'utiliser des moderations des fonctions de 
15 covariances. 

Uinteret d'utiliser ces modelisations est qu'elles permettent de 
disposer de matrices qui sont definies positives et inversibles. 

On comprend toutefois que ces techniques de filtrage, si elles 
donnent de bons resultats, sont fortement dependantes de I'expertise 
20 individuelle de la personne qui choisit les modelisations des differents 
variogrammes. 

II y a la une source d'erreurs, qui empeche que ces techniques 
puissent §tre utilis6es par des personnes qui n'en sont pas des specialistes. 

En outre, le choix des modelisations est egalement a Torigine de 
25 pertes de temps importantes en production. 

Un but de Tinvention est de pallier cet inconvenient et de proposer 
une technique de filtrage par analyse krigeante qui puisse etre mise en 
oeuvre de fagon automatique ou quasi-automatique. 

Elle propose un precede pour le filtrage d'au moins deux series de 
30 donnees sismiques representatives d'une meme zone, dans lequel on 
determine - par exemple par determination du variogramme crois6 de ces 
series de donnees et resolution de liquation de co-krigeage - une 
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estimation de la composante commune a ces series de donnees et on 
d§duit de cette estimation une decomposition de ces series de donnees. 

II est egalement propose un proc6d6 pour le traitement de donnees 
sismiques selon lequel - pour comparer deux series de donnees sismiques 
5 correspondant, pour une meme zone, a des grilles d'au moins un meme 
attribut obtenues pour deux valeurs distinctes d'au moins un param^tre 
donne - on met en oeuvre un proc6de de filtrage du type precite. 

Uinvention propose en outre un proc6d6 pour le filtrage d'au moins 
une s6rie de donnees representative des valeurs d'au moins un parametre 

10 physique sur au moins une zone, caracterise en ce qu'on identifie une 
mod6lisation d ! une composante de variabilite spatiale de son variogramme, 
on soustrait cette mod§lisation du variogramme experimental et on r6sout 
I'equation de krigeage correspondant a ce variogramme differentiel pour en 
deduire une estimation de la composante de variabilite correspondante sur 

15 la s6rie de donn§es. 

PRESENTATION DES FIGURES 

D'autres caracteristiques et avantages de Tinvention ressortiront 
20 encore de la description qui suit, laquelle est purement illustrative et non 
limitative et doit etre lue en regard des dessins annexes sur lesquels : 

- les figures 1 et 2a, 2b et 2c, dej& analysees, illustrent un exemple de 
cartographie de donnees sismiques et de decomposition par analyse 
krigeante correspondante ; 

25 - les 3a, 3b illustrent deux cartographies d'une meme zone, obtenues & 
partir d'acquisitions realis§es en deux temps differents ; 

- la figure 4 est une cartographie de la composante commune aux 
cartographie des figures 3a et 3b ; 

- les figures 5a et 5b et les figures 6a et 6b sont des cartographies des 
30 composantes autres que la composante commune des cartographies 

des figures 3a et 3b ; 

- les figures 7a a 7c sont des graphes illustrant la repartition des erreurs 
dans le cas respectivement d'un filtrage standard, d'un filtrage par 
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krigeage factoriel classique et d'un filtrage conforme a celui propose par 
Tinvention (co-krigeage factoriel). 

DESCRIPTION D'UN OU PLUSIEURS EXEMPLES DE MISE EN CEUVRE 
5 DE L'INVENTION 

Filtraae automatiaue 

On se place ici dans le cas oli on dispose de deux cartographies 
obtenues pour une meme zone avec des donnees sismiques acquises par 
10 exemple a des instants differents ou encore pour des attributs sismiques 
differents. 

Ces deux cartographies sont par exemple du type de celle illustr6e 
sur les figures 3a et 3b. 

Dans la suite du texte, on d6signera par Z1 et Z2 les deux fonctions 
1 5 correspondant a ces deux series de donnees. 

II est propose de decomposer chacune de ces deux fonctions en la 
somme de leur composante commune et de r6sidus orthogonaux. 

A cet effet, on determine a partir des deux series de donnees dont 
on dispose un variogramme crois6, dont les valeurs sont : 

20 Y 12 (h)= ~~r £ (Z1 (x) - Z1 (x+h))(Z2(x) - Z2(x+h)) 

N 

ou x et x+h designent les couples de points pris en consideration selon la 
direction et pour la distance h pour lesquels on determine la valeur du 
variogramme et 

oil N est le nombre de ces couples de points pour cette direction et cette 
25 distance. 

Connaissant ce variogramme crois6, on determine ensuite une 
estimation de la fonction qui lui correspond, laquelle v6rifie : 

30 oO a et p sont deux indices muets designant les points consideres autour du 
point x pour lequel on cherche a determiner une estimation de ladite 
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fonction, Z x a et Z 2 p etant la valeur audit point x, N etant le nombre de ces 
points et ou X x a et X 2 p sont des coefficients de ponderation. 

Ces coefficients de ponderation "k\ et \\ se determinent par 
inversion de i'equation de co-krigeage : 

5 
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ou les coefficients C12ap et C21aP sont les valeurs de variance 
croisee des fonctions Z1 et Z2 aux points correspondants aux indices a et 
pet ou les coefficients C11ap et C22ap sont les valeurs de covariance 
10 respectivement de la fonction Z1 et de la fonction Z2 auxdits points. 
.L'indice X correspond au point design^ par x prec6demment. 

On notera que la matrice qui apparaTt dans cette equation a 
Tavantage d'etre inversible sous certaines conditions de calcul. 

De cette faQon, en utilisant les covariances experimentales, on 
15 decompose de fapon automatique les deux variables correspondant aux 
deux series de donnees initiales en une composante commune et deux 
composantes residuelles orthogonales. La r§gularite des donnees fait que 
la covariance experimentale est connue pour toutes les distances utilis6es, 
aucune interpolation n'est n6cessaire, done la matrice est d6finie positive 
20 La fonction alors obtenue est une estimation de la composante 

commune aux deux series de donnees. 

La figure 4 represente une estimation de composante commune 
obtenue a partir des donnees correspondant aux cartographies des figures 
3a et 3b. 

25 On comprend qull est particulierement avantageux dans de 

nombreuses applications et notamment en sismique 4D de disposer de 
cette composante commune. 
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Elle permet notamment d'6liminer les artefacts d'acquisition non 
rep6tables 

- a partir de grilles d'attributs sismiques et notamment par exemple a 
partir d'amplitudes RMS dans un intervalle 

5 - a partir de grille de temps sismiques et par exemple a partir de temps 
d'evenements sismiques 

- a partir de volumes de vitesses sismiques et par exemple 3 partir de 
pointes automatiques de vitesses. 

Elle peut egalement etre utilisee en sismique simple pour eliminer 
10 des artefacts d'acquisition non r§p6tables 

- notamment a partir de grilles d'attributs sismiques calcules sur les 
classes d'incidence consecutives ou encore 

- £ partir de grilles d'attributs sismiques calcules sur des volumes issus de 
sommations partielles ou d'ondes converties. 

15 Par ailleurs, une fois cette estimation de la composante commune 

determinee, il est possible de determiner les composantes r§siduelles 
correspondant a la difference entre les donnees initiates et cette 
composante commune estimee, 

Ces composantes r6siduelles peuvent elles-memes etre 
20 decomposees par analyse krigeante. 

C'est ce qu'illustrent les figures 5a, 5b et 6a, 6b, qui sont des 
cartographies des bruits blancs et des bruits de raies lineaires estimes de 
cette fagon pour Tune et I'autre des deux series de mesures representees 
sur les figures 3a et 3b. 
25 Des exemples de mesures d f erreurs obtenues par filtrage standard, 

par filtrage par analyse krigeante puis par filtrage par analyse cokrigeante 
(ou krigejage multivariable) sont illustr6s sur les graphes des figures 7a & 7c. 

On comprend & la lecture de ces figures que le filtrage par analyse 
cokrigeante permet des dispersions bien moindres que dans les cas des 
30 filtrages classiques ou de filtrages par analyse krigeante et donne de 
meilleurs resultats. 

On a d6crit dans ce qui precede un exemple mis en ceuvre avec 
deux series de donnees. 
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Comme on le comprendra aisement, le procede propose peut 
egalement §tre mis en oeuvre de la meme fa?on avec un nombre de series 
de donn§es (campagnes) plus important. 

5 Filtraqe semi-automatique 

Ce deuxfeme mode de mise en oeuvre permet egalement une 
decomposition simplifiee dans le cas oil Ton ne dispose que d'une seule 
serie de donnees sismiques (fonction S1). 
10 II suppose que Ton dispose au prealable de la modelisation d'une 

composante em du variogramme experimental 9. 

Cette modelisation que Ton connaTt au pr6alable est par exemple la 
modelisation d'une composante independante de la geologie du sous-sol : 
bruit blanc, stries, etc. 
15 Connaissant cette modelisation de la composante 0m, on en deduit 

le variogramme residuel correspondant a la difference entre le variogramme 
experimental et cette composante 0m. 

On met alors en oeuvre une analyse krigeante pour determiner d'une 
part la composante modelisee Sm et d f autre part, a partir du variogramme 
20 r6siduel, le residu orthogonal R1 tel que : 
S1 = Sm + R1 

Cette decomposition automatique permet de filtrer les anomalies 
d'acquisition qui presentent une coherence spatiale ais6ment identifiable et 
modelisable, comme les stries parallfeles aux cables observees sur les 
25 amplitudes et les temps en sismique marine. 
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REVEND1CATIONS 

1. Precede pour le filtrage (fail moins deux series de donnees 
sismiques representatives d'une meme zone, caracterise en ce qu'on 

5 determine une estimation de la composante commune a ces series de 
donn6es et on d§duit de cette estimation une decomposition de ces series 
de donnees. 

2. Proc6d6 selon la revendication 1, caracterise en ce qu'on 
determine le variogramme croise de ces series de donnees et on resout 

10 Tequation de co-krigeage qui en resulte pour en deduire de fagon 
automatique une estimation de la composante commune a ces series de 
donnees. 

3. Precede selon Tune des revendication 1 ou 2, caracterise en ce 
que Ton determine des residus orthogonaux pour les differentes series de 

15 donn6es en soustrayant de chacune des series de donnees I'estimation de 
la composante commune. 

4. Precede selon Tune des revendications prec6dentes, caracterise 
en ce qu'on met en oeuvre une analyse krigeante pour decomposer lesdits 
r§sidus orthogonaux. 

20 5. Precede pour le traitement de donnees sismiques selon lequel, 

pour comparer deux series de donnees sismiques correspondant, pour une 
meme zone, a des grilles d'au moins un meme attribut obtenues pour deux 
valeurs distinctes d'au moins un parametre donne/on met en oeuvre un 
precede de filtrage selon Tune des revendications 1 a 4. 

25 6. Precede pour le filtrage d'au moins une s6rie de donnees 

sismiques representative d'au moins une zone, caracterise en ce qu'on 
identifie une modelisation d ! une composante de variability spatiale de son 
variogramme, on soustrait cette modelisation du variogramme experimental 
et on r£sout I'equation de krigeage correspondant a ce variogramme 

30 differentiel pour en deduire une estimation de la composante de variability 
correspondante sur la serie de donnees. 

7. Precede pour le traitement de donnees sismiques selon lequel, 
pour comparer deux series de donnees sismiques correspondant, pour une 
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meme zone, a des grilles d'au moins un meme attribut obtenues a deux 
instants differents, on met en ceuvre un proced§ de filtrage selon Tune des 
revendications 1 a 4. 
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